Exact fermion pair creation by an electric field in dS2 space-time 

S. Haouaijf] and R. Chekireb 

LPTh, Department of Physics, University of Jijel, 
BP 98, Ouled Aissa, Jijel 18000, Algeria. 

The effect of the electric field on the creation of spin 1/2 particles from vacuum 
in the (1+1) dimensional de-Sitter space-time is studied. The Dirac equation with 
a constant electric field is solved by introducing an unitary transformation. Then 
the canonical method based on Bogoliubov transformation is applied to calculate 
the pair creation probability and the density of created particles both for positive 
or negative wave vector. It is shown that in the presence of a constant electric field 
the density of created particles is more significant when the wave vector is in a well- 
defined direction. This depends on the orientation of the electric field and the sign of 
the particle charge. The limit H = 0, where dS space reduces to the flat Minkowski 
space-time, is studied. The imaginary part of the effective action is extracted from 
the vacuum to vacuum transition amplitude. 
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I. INTRODUCTION 

As we know our universe is undergoing an accelerated expansion that can be approxi- 
mated by an exponentially expanding de-Sitter space-time. It is widely believed, also, that 
the de-Sitter space-time allows a better description of the early stage of inflation. There- 
fore de-Sitter metric is of interest both for the cosmology of the early and late universe. 
This importance comes also from the fact that the de-Sitter space is the unique maximally 
symmetric curved space and enjoys the same degrees of symmetry as Minkowski space-time. 
This explains the increased interest to study physical quantum effects in such a space-time. 
Among these effects we cite the particle-antiparticle pair creation caused by the expansion of 
the universe [TH5]. The phenomenon of particle creation in de-Sitter space has been widely 
discussed and analyzed from various point of view [6HT5]. This is because the fact that the 
particle creation has many important applications in contemporary cosmology - e.g., it could 
have consequences for early universe cosmology and may play an important role in the exit 
from inflationary universe and in the cosmic evolution [Trj] - |2"T] . Also, the particle creation 
mechanism has significant contribution on the transition from an anisotropic universe to an 
isotropic one [22] . 

The purpose of this paper is to study the effect of the electric field on the creation of 
fermions from vacuum in the de-Sitter space-time. We know that actually there is no electric 
field in the universe. However, as is mentioned in [23], we reconcile that electric fields were 
present during the initial stages of the formation of the universe and they vanish because 
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of the inverse effect of the particle creation. This idea resembles to the assumption of the 
anisotropy of the early universe although the present universe is isotropic. Then the existence 
of such electric fields influences the cosmic evolution directly via Friedmann equations and 
by their effect on the creation of particles. 

We note that the influence of electric field on the particle creation in expanding universe 
has been studied in several scenarios describing different stages of the cosmic evolution 
[12 ] I22H28] . The creation of spin half fermions in (1+1) dimensional de-Sitter space-time by 
an electric field has been communicated several years ago in [29] where the author has found 
the positive and negative energy solutions for the corresponding Dirac equation. He has not 
shown however how those states are connected to one another by the charge conjugation. 
This symmetry between positive and negative solutions is indispensable in the derivation of 
the effect. 

In order to study the effect of particle-antiparticle pair creation form vacuum by grav- 
itational and electric fields we have at our disposal several methods such as the adiabatic 
method [301 IS], the Hamiltonian diagonalization technique [32] 133] . the Feynman path 
integral derivation [3U [35], the Green function approach [361 EZ], the semiclassical WKB 
approximation [38til0] as well as the method based on Bogoliubov transformation [HJ H2] 
that we shall use in this paper. The later method is more convenient for the spin 1 case. 

The paper is organized as follows; At the beginning we consider a spin 1 fermion sub- 
jected to an electric field in the (1 + 1) dimensional de-Sitter space-time. Then we solve 
the corresponding Dirac equation by introducing an unitary transformation. To investigate 
the process of particle creation we apply the canonical method based on Bogoliubov trans- 
formation connecting the " in" with the " out" states. This method permits us to determine 
the probability of pair creation, the density number of created particles and the vacuum 
persistence. To complete the analysis we write the probability of pair production per unit 
of time as a Schwinger series starting from the vacuum to vacuum transition amplitude and 
we show how the electric field amplifies particle creation. 



II. DIRAC EQUATION IN DE-SITTER SPACE 



Let us consider a spin | fermion of mass m and charge — e moves on the background 
geometry of the (1 + 1) dimensional de-Sitter space-time in the presence of an electric field. 
The line element of the metric describing dS2 space-time can be written as 

ds 2 = dt 2 - a 2 (t) dx 2 , (1) 

where a (t) = e Ht and H is the Hubble's constant. In terms the conformal time rj = 
J dt/a(t) = ^-e~ Ht this metric reads 

ds 2 = a 2 (77) {drf - dx 2 ) (2) 

where the new scale factor a (77) = a [t (r/)] = -=^. For the electric field we choose the gauge 
=( 0, A x (r/0), with 

This vector potential describes a constant electric field in the comoving system of coordi- 
nates. 
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The covariant Dirac equation with the metric ^ takes the form 

W (V) (fy - ieA^ (17) - I^fa) ) - m] V = (4) 
where the curvature-dependent Dirac matrices 7^(77) are given in diagonal tetrad gauge by 

where 7 and 7 1 are the usual Dirac's matrices that can be written, in (1+1) dimensional 
Minkowski space, in terms of Pauli matrices as follows 

7 = a z , 7 1 = ia y (7) 

and 1^(77) are the spin connections 

r = o 

2a(r/) 

By making the substitution x (Vi x ) = ^(rj)^ (77, x) , we obtain the simpler equation 

[7" (id^ + eA^rj)) - ma (77)] X (v, x ) = °- ( 9 ) 
In order to solve this equation we write, at the beginning, \ (77, x) = exp (ikx) £ (77) with 

where the two components £1 (77) and £ 2 (??) satisfy the two coupled equations 

Here we remark that the coupling coefficient k + depends on the conformal time and the 
usual iteration procedure leads to a complicated second order equation that does not admit 
well-known solutions. To simplify the problem let us introduce the unitary transformation 



6 (t?) J Vl + r 2 V ~ T 1 / V ^2 (rj) 



(13) 



with 



and 



T= 2L( M _ m ) (14) 
eEn 



M = \lmi+ e ^. (15) 
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Then the novel components ipi (rj) and if 2 (v) satisfy the following system of equations 

d 1 



drj 
. d 



MH 
1 



i^. 



eEnk 



eE k + 



M 2 Y 




m 


V ) . 


<P1 = 


M 


M 2 Y 




m 


V J. 


¥2 = 


M 



dr] MH 

which leads to a second-order equation for each component 



drj'- 



+ 



M 2 M 
lP ±l ll 



2eE n k 1 



n 



+ k 2 



kip 2 
k<fi 



0. 



(16) 
(17) 

(18) 



Let us notice that the gravitational field couples to the mass of the particle, while the electric 
field couples to the charge. In the novel system we can see that an effective field is coupled 
to the quantity M.. As we will show, for the process of particle creation, this quantity is 
more important than the mass of the particle. In effect, writing the equation (18) in the 



form (p" (rj) + u 2 (r])(p s (rj) = 0, with u 2 s (rj) = 
the adiabatic condition implies that 



^ + is — 

H 2 H 



1 + + k 2 , we can see that 



h 2 v 



lim 

r;-s>0 



U 2 S (V) 



H 
M 



« 1 



(19) 



Therefore the particle production is well-defined only when Ai » H. Furthermore, in the 
limit m — > 0, the mixing term in (16) and (17) vanishes and the positive and negative energy 
solutions never intercept each other. This means that there is no production of massless 
particles even if an electric field is present. 

In order to solve equation (18) we make the change of variable 77 — p, where 

p = —2iki]. 

The resulting equation is similar to the well-known Whittaker equation 



(20) 



Jp~ 2 



1 A 

i + P 



(p) = 0, 



(21) 



with (p s (p) = <p s (r)), s = 1,2 and the constants p 2 and A are given by 



/'•2 



1 M 

1 M 

2 +i H 



and 



A 



.eE 
'IP' 



(22) 



(23) 



It is known that one can find for the equation (21) several sets of linearly independent 
solutions which can be written in terms of the Whittaker functions M\ ifl (p) and W\ tll (p), 
with 



Mx,, (p) 
Wx,» (p) 



pH+?e~5M [ p 



p V+2 e -2TJ 



X+-,2p + l;p 



X+-,2p + l;p 



(24) 
(25) 
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where M (a, b, p) and U (a, b, p) are the Kummar functions jH]. 

The general technique for solving the Dirac equation in (1+1) dimensional de-Sitter 
space-time when an electric field is present being shown, let us use the obtained solutions 
to investigate the phenomenon of spin half particle creation. 



The study of particle creation in an expanding universe requires a definition of a vacuum 
state for the field theory [13]. However, unlike the free Dirac field, it is not obvious how to 
determine the vacuum states when the spinor field is subjected to a general gravitational 
background. It is widely believed that in arbitrary curved background, there is no absolute 
definition of the vacuum state and the concept of particles is not completely clear. From 
physical point of view it is well known that in the standard quantum theory a particle cannot 
be localized to a region smaller than its de Broglie wavelength. When this wavelength is 
sufficiently large, the concept particle becomes unclear [46J. Moreover, when the vacuum 
state is defined in the remote past it is habitually unstable so that it may differ from the 
vacuum state in the remote future. This gives rise to spontaneous particle creation. 

To obtain a well-defined vacuum state with a reasonable choice of positive and nega- 
tive frequency modes we use the so-called adiabatic method based on the solutions of the 
relativistic Hamilton-Jacobi equation. Taking into account the asymptotic behavior of the 
W\ tlx (p) functions 



and using the solutions of the Hamilton-Jacobi equation we can find for the " in" states the 
following components 



III. VACUUM DEFINITION AND PAIR CREATION 




(26) 




(27) 
(28) 



<fl,in (P) = (P) ' 



By the use of the relation |32] 




A W x ,,(p) = -(p + \ 



2 



1 



) 



Wx^-x (p) 



(29) 



we find the lower component for the positive and negative frequency modes 




(30) 



(31) 



For the Dirac spinor we have 



(V) = A/ - 



\-T^jM- e -fW- X ,,(-p) + 



■ yjM + fW-A,i-„ (-P) 



) 



(32) 
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and 



C„ (V) = A/" 



l M + (p) _ T M _ fW Ail _„ (p) 



Ty J M + (p ) + _ fW Ajl ^ (p) 

where M is a normalization constant that is unimportant vis-a-vis the mechanism of particle 



(33) 



creation. Note that we have multiplied {jj) bj^JM — ^jf- and £ in {jj) by y M + ^to as- 
sure that those states are connected to one another by the charge conjugation transformation 
defined by 

£ f = <xit- (34) 
In addition positive and negative energy solutions satisfy the orthogonality condition 



er = rt 



0. 



(35) 



With this choice we have a well-defined vacuum state which is no thing but the adiabatic 
vacuum. 

Let us now define the " out" states. As previously, by taking into account the asymptotic 
behavior of M A /J (p) at p — > 0, 



M x ^ (p) ~e-'/ + 2, 
we can define the following positive and negative frequency modes 

<fi,out (p) = m a,- m (p) 

<fI,out (p) = M -A,M (-P) • 

Then by the use of the functional relation 



(36) 



(37) 
(38) 



(2/i + l) 



g (2p + 1)' 
<9p 2p 



+ A 



^A, M (P) 



(P+^) 2 ~A 2 
2(p + l) (2p+l) 



we obtain 



in 

M 



out 



-M 



A-/H-1 



(p)- 



M A , M+1 (p) (39) 



(40) 



Consequently, we arrive at the following expressions for the two components of the Dirac 
spinor 



and 



Horn (v) = N* 



M 



A,— 



( P )+T 



M 



M 



A,-/u+l 



(P) 



-rM A ,_ M (p) + 



4 



M 



A,-/i+l 



(P) 



(41) 



(42) 



Using the fact that positive and negative energy solutions must be connected to one another 
by the charge conjugation we find that 



tlout (17) = ti h 



-P) + 



m 
M 



4(1 "if) 



M_A, t ("P) 



(43) 



7 



and 



&out (V) = A/" 



M_ 



-P) 



M 



4 



M_ 



(44) 



Now, as we have mentioned above, in order to determine the probability of pair creation 
and the density of created particles we use the Bogoliubov transformation connecting the 
"in" with the "out" states. The relation between those states can be obtained by the use 
of the relation between Whittaker functions H3l 



M x ,» (p) 



-iir\ 



r (2/x + 1) e 

r(/i-A + |) 



-w. 



-p) 



r (2/i + 1) e-M+D " 
+ w t r— n w ^ KP) 

r (/i + a + |j 



(45) 



with — | < argp < ^ and 2// 7^ —1, —2, • • •. The functions Cut (rj) and ^ t (rj) can be then 
expressed in terms of Cn (77) and (rj) as follows 



Cut (V) = <* Hn (V) + (3 £~„ (V) 

Cout (V) = «* Cn fo) + 0* fa) 

where the Bogoliubov coefficients a and /3 are given by 



(46) 



a A/" 



(47) 



and 



\a\' + \/3\ 2 = l. (48) 

The Bogoliubov relation between " in" and " out" states can be projected in the Fock space 
to be a relation between creation and annihilation operators 



b -kAn = «* h \out + P <out 



0>k,in — a &k,out + ft* b_ 



(49) 



k,ouf 



Therefore, the probability of pair creation and the density of created particles will be 
given in terms of Bogoliubov coefficients. In effect starting from the amplitude A = 

(0 out \a kjOUt b_ ktOUt \ in ), it is easy to show that 



A 



{®out\ Ojn) . 



The probability to create one pair of fermions in the state k is then 

2 



a 



By the use of the following property of gamma function 

12 



\r(ix) 



x sinh ttx 



(50) 



(51) 



(52) 
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we find 



sinh Ti ( ^ 



eEo\ 
H 2 I 



sinhvr^-f^ • (53) 
Notice that equation (53) differs from equation (45) in [29]. This is because the author of 
[29] used an erroneous formula between Whittaker functions (see equation (42) in [29J). 



For the density number of created particles we have 

n (k) = (0 in \al fiUt a kfiUt \ in ) 



(54) 



Taking into account that the Bogoliubov coefficients satisfy the condition (48), it is easy to 
show that 



n (k) 



sinh7r(f + ff)e" 



-2vr- 



(55) 



sinh7r (M _ <g>) +s inh7r (f + fa) e 
Let us note that this result is obtained for a positive wave vector {k > 0) and — f < argp < 



3- 



For the the case when k < 0, the quantities n (k) and Vk can be obtained form (53) 



and (55) by changing the sign of e. Since the particle creation is well-defined only in tne 



adiabatic limit, Ai >> H, the density of created particles can be approximated by 



n (k) = exp 



2vr 



M — sign (k) 



eEo 
H 



(56) 



which resembles to the Boltzmann distribution. This explains the thermal nature of the 



effect. In addition it follows from (56) that n (k) is more significant when k > 0. Therefore 
the constant electric field produces predominantly particles with k > 0. In other wards, in 
the presence of a constant electric field, particles prefer to be created in a specific direction. 
This depends on the orientation of the electric field and the sign of the particle charge. The 
antiparticles are in the main created in the opposite direction with the same density as the 
particles. 

At the end of this section let us examine the limit H — > where the de-Sitter space 
reduces to the flat Minkowski space-time. Taking into account that 4£ — ^ff ~ for 
small H, we easily obtain 



2eE ■ 



exp l-n (f ~ m „ 6XP H|) 
2sinh7r(f -ff) 



exp 



eE J 



and 



n (k) w exp I —n 



m 



eE C: 



(57) 



(58) 



We remark here that equations (57) and (58 ) are in complete agreements with the well-known 



results corresponding to the pair creation by an electric field in Minkowski space-time (see 
[4*8] and references therein). 



IV. SCHWINGER EFFECTIVE ACTION 

Let us now use the obtained results to calculate the imaginary part of the Schwinger 
effective action. We start by writing Vk as 
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where 



sinh7r(^ + e^)e - u 



sinh 7r ( ^ + e 
If 



#2 ; e 



exp 



M 



27T h + sinh 7r ( 



M 

H 



- ago ' 
" #2 , 



(60) 



where e = sign (fc). Let to be the probability to have no pair creation in the state k. The 
quantity CkVk is then the probability to have only one pair in the state k. Because of the 
Pauli principle we have Ck + CkVk = 1 and 

C k = 1 - a. (61) 

Next, we define the vacuum to vacuum transition amplitude by an intermediate effective 
action A (vac — vac) = exp (iS e ff). The vacuum to vacuum probability can be then written 
as | .A (vac — vac)\ = exp (— 2Im5 e //) = rife^- ^ follows from (61) that 



exp (—2 Im S e ff) = exp 



5>a 



and, consequently, 



2ImS, 



dk 



2na (t) 



By expanding the quantity ln(l — a), we obtain 



r 



dk 



2na (t) * — ' n 

v ; n=l 



E 1 



exp 



2vr 

H 



lnfl 



M 



' H 



(62) 



(63) 



(64) 



In the flat Minkowski space-time the integration over k can be carried out by making the 
change j dk — >■ cEq J dt = cEqT, which is consistent with the classical equation of motion 
cEq). In de-Sitter space, taking into account that the semiclassical approximation 



i.e. 



dk 
dt 



is valid as soon as the classical action is large, it follows from equation (A3) that the 



"in" and the "out" states behave like exponential functions with the classical action as 



exponent for >> 



eE$k 



and \t)\ « 



creation happens around the time r\ 



M 

eEok 
M 2 



respectively. This means that the particle 



, k or, equivalently, particles having a wave vector 
k are created on the average at the time rj ~ Al ~ 



eEgk ' 



set k 



M 



M z H_ e Ht 



Then, in de-Sitter space, we have to 

, L - ■ . Here t has to be understood as the time at which a particle is 

created with the wave vector k. Besides, at any time t, particles are produced with a wave 
vector in the range 

HM 2 m _ _ HM 2 



-k + 



<k<k + 



eE n 



(65) 



where kg is a positive constant. For the case when H ^ e Ht < kg we can write T in the form 



^ n 

n=l 

+E 1 



exp 



exp 



n=l 



2nn 

27171 



M 



ko- 



HM±Ht 



dk 



H 



M 



2na (t) 

dk 
ko+ HJ^ eHt 2na (t) 



exp 27m 



exp I —2-nn 



eEo 
H 2 

eE 



H 2 



(66) 
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The probability of particle creation per unit of time will be then given by 



dT _ m 2 H 2 + e 2 E 2 
dt eEn 



71=1 



27m 

IT 



M sinh 27m 



eEo 
H 2 



(67) 



This procedure allows us to write the number of created fermions per unit of time in the 
form 



dN _ m 2 H 2 + e 2 E 2 
dt neE 



sinh ( 27T 



eEo 
H 2 



cxp 



2tt 



M 



(68) 



which is the same as the first term in the Schwinger-like series (|67|) 

otice that 
that ^ and dr 



Notice that ^ and ^ are finite when the electric field is present. It is obvious also 



,,. <,.jim , t converge when E — >■ 0. The introduction of an electric field plays then 
an important role in the theory. This is a sort of regularization. When t is large we have 
H J£ e Ht > k and, consequently, the probability of particle creation per unit of time and 
the pair production rate become 

m 2 H 2 + e 2 E 2 ^ 1 

-E" ex P 



dr 

~dt 



eEn 



n=l 



n 



2im 



M 



eE 
H 



(69) 



and 



dN _ m 2 H 2 + e 2 E 2 
dt TteEo 



exp 



27m 



M 



eEp 
H 



As is done in the previous section when H — > equation (67) reduces to 



dT 

~dt 



eEo \ - 1 

27r ^ n 

n=l 



exp — nir 



m 
~e~E 



(70) 



(71) 



which is the same as equation (64) in 



V. CONCLUDING REMARKS 

In conclusion, we have studied in this paper the influence of an electric field on the 
creation of spin-^ particle pairs in the (1+1) inflationary de-Sitter space-time. We have 
found exact solutions for the corresponding Dirac equation. This is accomplished by the use 
of an unitary transformation which allows us to obtain a simple second-order differential 
equation for each component of the Dirac spinor. The charge symmetry between positive 
and negative frequency modes permitted us to find exact and analytic expressions for the 
Bogoliubov coefficients without need to normalize the wave functions. Then the pair creation 
probability and the density of created fermions are calculated by applying the canonical 
method based on Bogoliubov transformation connecting the "in" with the "out" states. 
From these results we were able to extract the imaginary part of the Schwinger effective 
action. 



We have found that equation ( 53 ) in the present paper differs from equation (45) in 
This is due to a mistake in equation (42) in [29J which leads to incorrect results. Furthermore, 
equation (46) in [29J is obtained by considering the case when \ § ± ^| » 1 and using the 
approximation sinli7r (7^+ 77?) ~ sinli7r (7^ — 7^) which is valid only for a weak electric 
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field. When the effect of the electric field is significant (i.e. the case of strong field) the good 
approximation can be written as sinli7r (j£ + ^) « exp (27r^) sinh7r (~ — Tp 1 )- 

We conclude through this study that a constant electric field amplifies the creation of 
fermions with positive wave number k and minimizes it in the opposite direction. The inverse 
is true for antiparticles. In addition we have shown that the creation of light particles is 
possible in the presence of the electric field and thermal spectrum of particle creation will be 
obtained only when M. satisfy the adiabatic condition j£ >> 1. Another important result 
is that the creation of massless particles with conformal coupling is impossible even if an 
electric field is present. 

Let us note at the end of this work that the method used to solve the Dirac equation 
is useful not only in the present case but also for all situations where A z (77) has the same 
dependence on 77 as a (77). This encourages us to consider more realistic cases and especially 
the (3+1) dimensional problem. 



Appendix A: The relativistic Hamilton-Jacobi equation 

The " in" and " out" states are chosen according to the asymptotic behavior of the semi- 
classical solutions constructed by the use of WKB approximation starting from the solutions 
of the relativistic Hamilton-Jacobi (H-J) equation. It is widely known that in the adiabatic 
condition, i.e. M. » H, this approximation proved most fruitful in studying the particle 
creation either for scalar particles or for fermions. In this appendix we give the analytic 
solution of the H-J equation. The general form of the H-J equation is given by 

<T (d^S + eAJ (d u S + eA u ) - m 2 = (Al) 

For the case of dS2 space with constant electric field the classical action S can be decomposed 
as 

S = G (77) + kx, (A2) 
where the time dependent part G (77) satisfies the following equation 



By the use of the well-known integral (13 

y/a + fix + 7X 



2 

dx 



x 



I 7, n 1- ■ ! f 2a + f3x 

\? a + px + *yx z — y/a arcsmh — = 

where a, (3 and 7 are real numbers, with 7 > and 4a7 — (3 2 > 0, we obtain 



12 




Notice that at rj — > 0, we have G (rj) ±.Mt+a Constant. 
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